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Tropical optimization problems in 
time-constrained project scheduling* 

N. Krivuliit 


Abstract 

We consider a project that consists of activities to be performed in 
parallel under various temporal constraints, which include start-start, 
start-finish and finish-start precedence relationships, release times, dead¬ 
lines, and due dates. Scheduling problems are formulated to find opti¬ 
mal schedules for the project with respect to different objective func¬ 
tions to be minimized, such as the project makespan, the maximum 
deviation from the due dates, the maximum flow-time, and the max¬ 
imum deviation of finish times. We represent these problems as op¬ 
timization problems in terms of tropical mathematics, and then solve 
them by applying direct solution methods of tropical optimization. As 
a result, new direct solutions of the scheduling problems are obtained 
in a compact vector form, which is ready for further analysis and prac¬ 
tical implementation. The solutions are illustrated by simple numerical 
examples. 

Key-Words: idempotent semifield, optimization problem, project 
scheduling, precedence relationship, scheduling objective. 
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1 Introduction 

Tropical optimization problems, which are formulated and solved in the 
framework of tropical mathematics, find increasing use in various fields of 
operations research, including project scheduling. As an applied mathe¬ 
matical discipline concentrated on the theory and applications of semirings 
with idempotent addition, tropical mathematics dates back to a few semi¬ 
nal papers [331 El da Elsa S3 ED], including those mm concerned with 
optimization problems drawn from machine scheduling. 
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In succeeding years, tropical mathematics was studied by many authors 
under various names, such as idempotent algebra, max-algebra, max-plus 
algebra, extremal algebra (see m usi im [i6j m [ 31 ] esi Esi m and references 
therein). Tropical optimization problems were investigated in a number of 
works, in which scheduling issues frequently helped to motivate and illustrate 
the study. Specifically, tropical mathematics (max-plus algebra) served as a 
solution framework for scheduling problems in mmmmmummm 
EU32U2]. 

Many optimization problems are formulated in the tropical mathematics 
setting to minimize nonlinear functions defined on vectors over idempotent 
semifields (semirings with multiplicative inverses), subject to constraints 
given by vector equations and inequalities (see, e.g., overviews in mm)- 
For some problems, direct solutions are obtained in a closed form under 
general assumptions. Other problems are solved algorithmically by using 
iterative computational procedures. 

This paper examines problems that are drawn from time-constrained 
project scheduling, which involves the planning of activities in a project 
over time to achieve certain objectives PU E3 SUES]. The aim of the study 
is to provide new solutions to the problems on the basis of recent results in 
tropical optimization. 

We consider scheduling problems, which are to find an optimal schedule 
for a project that consists of a set of activities operating in parallel under 
various temporal constraints, including start-start, start-finish, finish-start, 
release time, deadline, and due-date constraints. As optimization criteria 
to minimize, we take the project makespan, the maximum deviation from 
due dates, the maximum flow-time, and the maximum deviation of finish 
times. Such problems are known to have algorithmic solutions in the form 
of iterative computational procedures. Specifically, many problems can be 
formulated and solved as linear, integer or mixed-integer linear programs 
by appropriate mathematical programming algorithms. Examples of the 
algorithmic solutions can be found in [MUSE] (see also reviews in P3JE2 
SUES]). 

We represent the scheduling problems as tropical optimization problems, 
which are then solved by applying solution methods developed in [2j 2? 26]. 
We derive new direct solutions to the scheduling problems considered, which, 
in contrast to the conventional algorithmic solutions, provide results in a 
compact explicit vector form, ready for further analysis and applications, 
and thus have the potential to complement and supplement existing ap¬ 
proaches. These solutions allow various constraints to be incorporated in a 
unified and constructive way. The calculation of the solutions involves sim¬ 
ple matrix and vector computations according to explicit formulae, which 
offers a basis for the development of efficient computational algorithms and 
software implementation. 

The current paper further extends and improve the results presented in 
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the conference paper [28] in an effort to incorporate illustrative numerical 
examples, discuss computational complexity of solutions, enhance formulae, 
and refine the reference list. In the paper, we continue the research on 
the application of tropical optimization to scheduling problems reported 
in journal and conference publications mmmm [29] . We examine 
more complicated general scheduling problems with new objective functions 
introduced and additional temporal constraints imposed. We show how 
tropical optimization techniques can be applied to these problems, which 
results in new solutions in both tropical optimization and project scheduling. 

The rest of the paper is organized as follows. In Section [D we describe 
scheduling problems that motivate and illustrate the study, and formulate 
these problems by using the conventional notation. Section [3] includes a brief 
overview of preliminary definitions and results of tropical mathematics to 
be used in the subsequent sections. In Section |4j we present some tropical 
optimization problems and their solutions, and discuss the computational 
complexity of the solutions. In Section [5] we first rewrite the scheduling 
problems as tropical optimization problems, and then solve them by applying 
the results from Section [3J The solutions are illustrated by simple, but 
representative, numerical examples. 

2 Project scheduling model and example problems 

We start with the description of a project scheduling model in a general 
form, and then present example problems to find optimal schedules. To 
facilitate the subsequent representation of the problems in terms of tropical 
mathematics, we employ a somewhat different notation than that commonly 
adopted in the literature (see, e.g., hd m\ eii m for further details and 
standard notations of project scheduling). 

Consider a project that consists of n activities operating in parallel under 
start-start, start-finish and finish-start precedence relations, due dates, and 
time boundaries for start and finish times in the form of release time, release 
deadline and deadline constraints. To describe the temporal constraints and 
scheduling objectives under consideration, we use the symbols Xi and y*, 
which, respectively, represent the unknown start and finish times for each 
activity i = 1 ,... , n. 

2.1 Temporal constraints 

We now examine constraints imposed on the start and finish times of each 
activity i = 1, ...,n. First, we represent precedence relations, which link 
activity i with other activities. Let a^- be the minimum possible time lag 
between the start of activity j and the finish of i . The time lag an specifies 
the minimum duration of the activity (the duration provided that no other 
constraints are imposed). Note that the value — a^- can be interpreted as 
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the maximum time lag between the finish of j and the start of i. If there 
is no lag defined, we assume = — oo. 

The start-finish constraints take the form of the inequalities + Xj < yi 
holding for all j = 1 ,,n. The activity is assumed to finish immediately 
after all related start-finish constraints are satisfied, and thus at least one of 
the inequalities must hold as an equality. Then, these inequalities combine 
to give one inequality, which, under the above assumption of immediate 
finish, is equivalent to the equality 

max (aij + Xj) = yi. 

l<jf<n 

Furthermore, we use the notation bjj to describe the minimum time lag 
between the start of activity j and the start of i , and put bij = —oo if the 
lag is not specified. The start-start constraints are given by the inequalities 
bij + Xj < Xi for all j , which can readily be rewritten as one inequality 

max (bn + Xj) < Xj. 

1 <j<n J J 

Let the minimum time lag between the finish of activity j and the start of 
i be denoted by Cjj , with Cjj = —oo if undefined. The finish-start constraints 
are written as the inequalities Cjj + yj < Xi for all j , or as one inequality 

max (dj + yj) < Xj. 

1 <j<n J J 

Finally, we introduce due dates and time boundary constraints. The due 
date indicates the time when the activity is ideally expected to finish. Since 
the due date may be unachievable under other constraints, it is considered 
as not a strict constraint. For activity i , we denote the due date by dj. 

Let gj and hi be the earliest and latest possible times to start, and /j 
be the latest possible time to finish. The release time, release deadline, and 
deadline constraints provide strict lower and upper boundaries for the start 
and finish times, given by 

9i<Xi< hi, yi < /j. 

2.2 Optimization criteria 

To describe scheduling objectives, we use several criteria that commonly 
arise in the development of optimal schedules in practice. The criteria are 
written below in the form, which is ready for immediate translation into 
terms of tropical mathematics. 

We begin with the maximum absolute deviation of finish times of ac¬ 
tivities from due dates that a project should meet. The minimum value of 
this criterion corresponds to the least violation of the due dates, which can 
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be attained. With the notation introduced above, the maximum deviation 
from the due dates is given by 

max |yi — di\ = max max(;% -di,di~ m). 

1 <i<n l<i<n 

Next, we consider the maximum deviation of completion times of all ac¬ 
tivities. The minimization of this criterion is equivalent to finding a schedule, 
where all activities have to finish simultaneously as much as possible. Such 
a problem can arise in just-in-time manufacturing, when certain delivery op¬ 
erations must be completed at once. The maximum deviation of completion 
times is written as 


max yi — min % = max % + max (—%). 
l<i<n l<i<n 1 <i<n l<i<n 

The flow-time of an activity (also known as the system, throughput and 
turn-around time) is defined as the difference between its start and finish 
times, and can determine expenses related to undertaking the activity in a 
project. The flow-time of activity i is bounded from below by the value 
of an , which is commonly assumed to be nonnegative, and may be greater 
than an due to other temporal constraints. 

In many real-world problems, the objective is formulated to minimize 
the maximum flow-time taken over all activities, and thus described by the 
expression 

max (yi — Xi). 
i<*<« 

Finally, we discuss the makespan, which is the interval between the ear¬ 
liest start time and the latest finish time of activities in a project. The 
makespan indicates the total duration of the project, and finds wide applica¬ 
tion as an objective function to be minimized in many scheduling problems. 
The makespan is given by 


max yj — min a%- = max in + max (—xs). 
l<i<n Ki<n Ki<n \<i<n ’ 


2.3 Examples of scheduling problems 

We conclude with typical examples of scheduling problems, which are to 
serve to both motivate and illustrate the results in the rest of the paper. 
To formulate the problems, we use the notation and formulae introduced 
above for the unknown variables, given parameters, temporal constraints 
and scheduling objectives. 

2.3.1 Minimization of maximum deviation from due dates 

First, we consider a problem to minimize the maximum deviation from due 
dates under start-finish, start-start and finish-start constraints. Given the 
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parameters a t j , b tJ , Cy and d, , the problem is to find the unknown start 
time ay and finish time y* for each activity i = 1,... ,n, that 


minimize 

max max(t/ ? - — 

l<j<n 

di , di 

- Vi), 


subject to 

max (an + ay) 
i<j<n 

= yu 

max (bn + ay) < ay, 
l<j<n 

(1) 


max (cy + yj) 
l<7<n 

< Xi, 

i = 1 ,..., n. 



2.3.2 Minimization of maximum deviation of finish times 

We now formulate a problem of minimizing the maximum deviation of fin¬ 
ish times, subject to start-finish, start-start, finish-start, and deadline con¬ 
straints. Given the parameters «y, bij, Cy and fi, we need to determine 
the unknowns ay and y* that 


minimize max y,• + max (—y, 

1 <i<n 1 <i<n ' 

subject to max (oy + Xj) = y,;, max (b tj + Xj) < Xi, 
l<j<n l<j<n 

max (cij + y,) < ay, yi < fi, * = 1 ,..., n. 

l<j<n 


( 2 ) 


2.3.3 Minimization of maximum flow-time 

Next, we consider a problem of minimizing the maximum flow-time under 
start-finish, start-start, finish-start and release time constraints. Given the 
parameters a VJ , 6y, c tJ and y*, we find the values of Xi and y* that solve 
the problem 

minimize max (y, — ay), 
l<i<n 

subject to max (an + ay) = y,;, 

l<j<n 

max (c™ + yj) < ay, 

l<j<n 


max (bin + Xn) < ay 
1<j<n J J 


(3) 


y* ^ J'i, 


* = l, 


,n. 


2.3.4 Minimization of makespan 

Suppose that we need to minimize the makespan of a project subject to start- 
finish, release time, release deadline, and deadline temporal constraints. 
Given ay, yj, hi and fi , we find ay and yi that 


minimize max yi + max (—ay), 

l<i<n l<i<n 

subject to max (a,;, + ay) = y*, 

l<j<n 

gi<Xi< hi, yi< fi, i = l,...,n. 


(4) 
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Note that problems like those presented above can normally be solved us¬ 
ing iterative computational algorithms (see, e.g., UDESIETIES] for overviews 
of available solutions). Specifically, these problems can be formulated as 
linear programs to solve them by computational methods of linear program¬ 
ming, which generally offer algorithmic solutions. Below, we provide new 
solutions to the problems, which are based on optimization methods in tropi¬ 
cal mathematics, and present results in a compact explicit vector form rather 
than in the form of a numerical algorithm. 

3 Preliminary algebraic definitions and results 

In this section, we give a brief overview of preliminary definitions, notation 
and results of tropical algebra to provide a formal basis for the description 
and application of tropical optimization problems in the next sections. Both 
introductory and advanced material on tropical mathematics can be found 
in many publications, including Bi nan use Buna Isa ei to name only a 
few. The overview given below is mainly based on the presentation of results 
in [2H [271 [26], which offers a useful framework to obtain direct solutions to 
the problems under study in a compact vector form. 

3.1 Idempotent semifield 

Consider a system (X, ©, ©, 0,1), where X is a set, which is closed under 
addition © and multiplication © with zero 0 and identity 1, such that 
(X, ©, 0) is a commutative idempotent monoid, (X \ {©}, ©, 1) is an Abelian 
group, multiplication is distributive over addition, and © is absorbing for 
multiplication. This system is usually called the idempotent semifield. 

Addition is idempotent, which means that x © x = x for each x £ X. 
The idempotent addition induces on X a partial order such that x < y if and 
only if x © y = y. It follows directly from the definition that x < x © y and 
y < x®y for all x, y € X. Moreover, the inequality x®y < z appears to be 
equivalent to the two inequalities x < z and y < z, and both addition and 
multiplication are monotone in each argument. Finally, the partial order 
extends to a linear order on X. 

Multiplication is invertible to let each nonzero x £ X have the inverse 
x _1 such that x©x _1 = 1. The inverse operation is antitone, which implies 
that, for all nonzero x and y , the inequality x < y yields x _1 > y _1 . 

The power notation with integer exponents is routinely used to represent 
repeated multiplication, and defined as follows: x° = 1, x p = x © x p_1 
and x~ p = (x -1 ) p for all nonzero x and positive integer p. The integer 
powers are assumed to extend to rational exponents to make X algebraically 
complete. 

In the algebraic expressions below, the multiplication sign © is omited 
to save writing, and the exponents are read in the sense of tropical algebra. 
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An example of the idempotent semifield under consideration is the real 
semifield R m ax,+ = (RU {— oo}, max, +, — oo, 0), in which the addition © is 
defined as maximum, and the multiplication © is as ordinary addition, with 
the zero 0 given by —oo, and the identity 1 by 0. Each number x G 1R has 
the inverse x _1 equal to the opposite number — x in the conventional nota¬ 
tion. For all x,y G R, the power x y is well-defined and coincides with the 
arithmetic product xy . The partial order induced by idempotent addition 
corresponds to the standard linear order on [R. 

3.2 Matrices and vectors 

We now examine matrices and vectors over the idempotent semifield intro¬ 
duced above. The set of matrices that have m rows and n columns with 
entries from X is denoted X mxn . A matrix with all entries equal to 0 is the 
zero matrix. A matrix is called row-regular (column-regular), if it has no 
rows (columns) that consist entirely of 0. Provided that a matrix is both 
row- and column-regular, it is regular. 

For any matrices A, B G X mxn and C G X nx ^, and a scalar x G X, 
the matrix addition, matrix multiplication and scalar multiplication follow 
the standard rules with the scalar operations © and © in the place of the 
ordinary addition and multiplication, and are given by the formulae 

n 

{A © B}ij = {A}ij © {B}ij, {AC}ij = (^{A}ik{C}kj, 

k= 1 

{x A ) ,y — x{A}ij. 

The partial order associated with the idempotent addition and its prop¬ 
erties extend to the matrices, where the relations are expanded entry-wise. 

For any matrix A = ( a\j ) G X mxn , the transpose is the matrix A T G 

y^nxm 

The multiplicative conjugate transpose of A is the matrix A~ = (a“) G 
X nxm with the entries a© = a© 1 if a,-; © 0, and a© = 0 otherwise. 

J L J 

Consider square matrices of order n in the set X nxn . A matrix having 
the diagonal entries equal to 1, and the off-diagonal entries to 0, is the 
identity matrix denoted I. The power notation with nonnegative integer 
exponents serves to represent iterated products as follows: A 0 = I and 
A p = A P ~ 1 A for any matrix A and integer p > 0. 

For any matrix A = (a,y) G X nxn , the trace is given by 

n 

tr A — cqi © * * * © Oj nn — ^J^) eta- 

i =1 

For any matrices A and B, and scalar x, the following identities are 
valid: 


tr(A © B) = tr A © tr B, tr(AB) = tr(BA), tr(xA) = xtr A. 



Every matrix having only one row (column) is considered a row (column) 
vector. All vectors below are column vectors unless otherwise specified. The 
column vectors of order n form the set X n . A vector with all elements equal 
to 0 is the zero vector. If a vector has no zero elements, it is regular. The 
vector of all ones is 1 = (1, ..., t) T . 

Let A £ X nxn be a row-regular matrix and x £ X n a regular vector. 
Then, the vector Ax is regular. If A is column-regular, then the row vector 
x T A is regular. 

The multiplicative conjugate transpose of a nonzero column vector x = 
(xi) G X" is a row vector x~ = (x“) with the entries = a;” 1 if Xi / 0, 
and Xi = 0 otherwise. 

The conjugate transposition has the following useful properties. First, 
for any nonzero vector x, the equality x~x = 1 is valid. 

Furthermore, if x and y are regular vectors of the same order, then 
the element-wise inequality x < y is equivalent to x~ > y~. In addition, 
the matrix inequality xy~ > (x~y)~ l I holds, and becomes the inequality 
xx~ > I when y = x. 

Finally, consider a square matrix A £ X nxn . A scalar A £ X is an 
eigenvalue of A, if there exists a nonzero vector x £ X n to satisfy the 
equality Ax = Xx. The maximum eigenvalue of the matrix A is called the 
spectral radius, and given by 

n 

A = tr A © • • • ® tr V n (A n ) = 0 ti l/k (A k ). 

k =1 


3.3 Solution to linear inequalities 

We conclude the overview of the preliminary results with the solutions to 
vector inequalities, which arise in the literature in different settings and 
have solutions given by many authors in various forms (see, e.g., mm- 
Below, we describe explicit solutions represented in a compact closed form 
that provides a unified framework for the systematic analysis of optimization 
problems in what follows. 

Suppose that, given a matrix A £ X mXTl and a regular vector d £ X m , 
we find vectors x £ X n to solve the inequality 

Ax < d. (5) 

A direct solution, which uses algebraic properties of the multiplicative 
conjugate transposition to reduce © to an equivalent inequality, is obtained 
in [26] as follows. 

Lemma 1. For any column-regular matrix A and regular vector d, all 
solutions to inequality © are given by 

x < (dT A)~. 
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Furthermore, we consider the problem: given a matrix A € X nxn , find 
regular vectors x £ X n that satisfy the inequality 

Ax < x. (6) 

To represent a solution to the problem for any matrix A £ X nxn , we 
define a function that takes A to the scalar 

n 

Tr(A) = tr A © • • • © tr A n = ®trA fc , 

fc=i 

and make use of the asterisk operator (also known as the Kleene star), which 
maps A with Tr(A) < 1 to the matrix 

n— 1 

A* = I © A © ■ ■ ■ © A”" 1 = 0 A k . 

k =0 

The next result, based on the properties of the asterisk operator, is 
derived in [26] to offer a complete solution to inequality © (see also ©HE?]). 

Theorem 2. For any matrix A, the following statements hold: 

1. If Tr(A) < t, then all regular solutions to © are given by x = A*u, 
where u is any regular vector. 

2. If Tr(A) > t, then there is no regular solution. 

The results, given by Lemma [1] and Theorem [2] provide necessary in¬ 
struments for the solutions of optimization problems presented in the next 
section, as well as for the application of the solutions to scheduling problems 
in the last section. 

4 Tropical optimization problems 

Tropical optimization problems present an area in tropical mathematics, 
which is of both theoretical interest and practical importance (see, e.g., 
[S S2l El |25l for further details and overviews). Many problems are for¬ 
mulated in the framework of tropical mathematics to minimize nonlinear 
functions defined on vectors over idempotent semifields and calculated using 
multiplicative conjugate transposition of vectors. These problems may have 
constraints, which are given by vector equations and inequalities. There are 
problems that can be solved directly in a general setting. For other problems, 
only algorithmic solutions are known, which offer an iterative computational 
scheme to produce a solution, or signify that no solutions exist. 

The purpose of this section is twofold: first, to offer representative ex¬ 
amples to demonstrate a variety of optimization problems under study, and 
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second, to provide an efficient basis for the solution of scheduling problems 
in the next section. We consider examples of both unconstrained and con¬ 
strained optimization problems with different objective functions defined in 
the common setting in terms of a general idempotent semiheld. For all 
problems, direct solutions are given in a compact vector form ready for fur¬ 
ther analysis and straightforward computations. For some problems, the 
solutions obtained are complete solutions. 

4.1 Examples of optimization problems 

We begin with the following problem. Suppose that, given a matrix A £ 
X mxn and a vector d £ X m , we need to find regular vectors x £ X n that 

minimize d~ Ax © (Ax)~d. (7) 

The next statement offers a direct algebraic solution to the problem, 
which is obtained by deriving a strict lower bound for the objective function 
together with a particular solution that yields the bound (see, e.g., [22]). 

Theorem 3. Let A be a row-regular matrix and d a regular vector. Then, 
the minimum value in problem © is equal to 

A = (( A(d-A)-)~d) 1 ' 2 , 

and the maximum solution is given by 

x = A (d~A)~. 

Furthermore, suppose that, given matrices A,B £ X mxn and vectors 
p,q £ X m , the purpose is to obtain regular vectors x £ X n to solve the 
problem 

minimize q~Bx(Ax)~p. (8) 

A direct solution to the problem can be found, using a similar technique 
of evaluating a strict lower bound as above, in the following form |23j. 

Theorem 4. Let A be row-regular and B column-regular matrices, p be 
nonzero and q regular vectors. Then, the minimum value in problem ([8j) is 
equal to 

A = ( A(q~B)~)~p , 

and attained at any vector 

x = a(q~ B)~ , a > 0. 


11 


Given matrices £ X" x " and a vector g G X n , consider the problem 

to find regular vectors x G X n that 

minimize x~ Ax , 

(9) 

subject to .Bai © g < s. 

A direct solution to the problem is given in [241 [27] by introducing an 
auxiliary parameter to represent the minimum value of the objective func¬ 
tion, and reducing the problem to a parameterized vector inequality. The 
existence condition for solutions of the inequality is used to evaluate the 
parameter, whereas the solutions of the inequality are taken as the complete 
solution to the optimization problem. 

Theorem 5. Let A be a matrix with spectral radius A > 0 , and B be a 
matrix such that Tr (B) < 1. Then, the minimum value in problem Q is 
equal to 

n— 1 

0 = A © 0 0 tv 1/k ( AB h ■■■AB ik ), 

k= 1 l<iiH- \-ik<n—k 

and all regular solutions are given by 

x = ( 0” 1 A © B)*u , u > g. 

Finally, suppose that, given a matrix A G X nxn and vectors g,h G X n , 
we need to find regular vectors that solve the problem 


minimize x Ax, 
subject to g < x < h. 


( 10 ) 


The same technique as above yields the next solution proposed in [24] . 

Theorem 6. Let A be a matrix with spectral radius A > 0, and h be a 
regular vector such that h~g < t. Then, the minimum value in problem 
m is equal to 

n— 1 

Q = \ ®(^(h~A k g) l / k , 

k=1 

and all regular solutions are given by 

x = (8~ 1 A)*u, g < u < (h~ (O^ 1 A)*)~. 


4.2 Computational complexity of solutions 

We conclude this section with some remarks on the computational complex¬ 
ity of the results with respect to the dimension n of the problems. First, 
note that the computation of the solutions provided by Theorems [3] [7] and [6] 
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to problems 0, © and m involves only a fixed number of basic matrix- 
vector operations, and thus obviously has a polynomial complexity. Specif¬ 
ically, the calculation of the star A* for a matrix A of order n requires 
computing of the sum of the n — 1 first powers of A , and therefore, takes no 
more than 0(n 4 ) scalar operations. The computation of the spectral radius 
A for the matrix has the same level of complexity. 

Consider the solution given by Theorem 0 to problem ([U]) , and verify 
that it can be calculated in polynomial time. We write the minimum in the 
problem as 

n— 1 n—k 

0 = A©//, /x = ^^tr^^Tfcz), 

k =1 1=1 

where TJy is the sum of matrix products AB n ■ ■ ■ AB lk over all nonnegative 
integers ii ,..., ik such that i\ + • • • + ik = l. 

The matrices Tjy satisfy the recurrence equation T^j = Tk~\j,A®Tkj,-\B , 
where Tk o = A k , Tqi = B 1 and Too = T, which involves one matrix addition 
and two matrix multiplications per matrix. Since the number of matrices 
Tf-i needed to evaluate the sum /x is 1 + ■ ■ ■ + n = n(n + l)/2, the complex¬ 
ity of computing /x, and hence of 9, is at most 0(n 5 ). Moreover, given the 
value of 9, any solution vector x is computed in polynomial time as well, 
and thus we conclude that the overall solution has polynomial complexity. 

5 Application to project scheduling 

We are now in a position to derive solutions to the scheduling problems 
formulated in the beginning of the paper. In this section, we first represent 
each problem in the framework of the idempotent semifield lR m ax,+ in both 
scalar and vector forms, and then solve it by reducing to an optimization 
problem of the previous section. 

To provide insight into the solution method proposed, and to illustrate 
the computation technique used, we offer numerical examples of solving 
scheduling problems. To save room, we restrict the examples to artificial 
problems of low dimension, which, however, unambiguously demonstrate the 
practicability of the approach to solve real-world problems of high dimension. 

5.1 Minimization of maximum deviation from due dates 

Consider problem (|T]) and describe it in terms of the semifield tR max ,+ • By 
replacing the usual operations by those of lR ma x,+ , we obtain the problem 
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to find the unknowns Xi and r/i for all i = 1,... ,n to 

n 

®K _1 y* ®y~ l di), 

i =1 

n n n 

® aijXj = yi, ® bijXj < Xi, ® Cijyj < x i: 

3 = 1 i =1 i =1 

i = 1,... ,n. 

To put the problem in a vector form, we introduce the matrix-vector 
notation 

A = ( a i:j ), B = ( bij), C = ( Cij), d = ( d t ), x = (xi), y = ( yi ). 

With this notation, the problem is to find the unknown vectors x and 
y that 

minimize d~y®y~d, 

subject to A.x = y, Bx ^ x, Cy < X. 

The next result offers a solution to the problem. 

Theorem 7. Let A be a row-regular matrix, the matrix D = B © CA 
satisfy the condition Tr(T>) < 1, and d be a regular vector. Then, the 
minimum value in problem cu is equal to 

A = ((AD*{d~AD*)~yd) 1/2 , 

and the maximum solution is given by 

x = AD*{d~AD*)~, y = AAD*(d~AD*)~. 

Proof. After substitution y = Ax, we combine both inequality constraints 
into one inequality of the form Dx < x with D = B®CA. An application 
of Theorem [2] to this inequality yields the solution x = D*u, where u is 
any regular vector. 

Substitution of this solution reduces problem (fTTj) to the unconstrained 
problem 

minimize d~ AD*u © (AD*u)~d. 

The last problem has the form of © with A replaced by AD*. There¬ 
fore, we can apply Theorem [3] to obtain a solution in terms of the unknown 
vector u. Turning back to the vectors x and y leads to the desired re¬ 
sult. □ 

We now consider the conditions of the theorem to see that, in the con¬ 
text of project scheduling, these and similar conditions are naturally met in 
practice, unless the problem under consideration is incorrectly posed, due 


minimize 


subject to 
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to wrong or incompatible conditions. Specifically, since the elements on 
the diagonal of the matrix A represent the minimum duration of activities, 
and thus must be greater than 0 = —oo, the matrix is formally both row- 
and column-regular. Moreover, for the same reason, the matrix A has the 
spectral radius A > 0. 

The regularity assumption on the vector d is equivalent to assuming 
that the deadlines for all activities are finite (not equal to 0 = —oo), which 
is the case for the real-world problems, and thus this assumption is typically 
fulfilled. The same conclusion can be reached for similar vectors encountered 
in the subsequent proofs. 

Finally, the condition Tr(IA) < 1 implies that the constraints in the 
problem are compatible to provide nonempty feasible sets for the unknown 
vectors x and y. 

We illustrate the direct solution offered by Theorem [7] by the following 
example. 

Example 1. Consider a project that involves n = 3 activities that operate 
under start-finish, start-start, finish-start and due dates constraints, given 
by the following matrices and the vector: 



where the symbol 0 = —oo is used to simplify the writing of matrices. 

To verify the conditions of Theorem [71 we first note that the matrix A 
is row-regular and the vector d is regular. Next, we calculate the matrices 

/ - 1 - 3 1 \ 

CA = 1-13, D = B ®CA = 

\ 0 0 0 / 

and then take the matrix D to find the powers 

/ 0 -3 1 \ / 0 -2 1 \ 

D 2 = 2-12, D 3 = I 1 03. 

\ -2 -3 0 / \ —1 —4 0 / 

After evaluating the traces, we have Tr (D) = tr D ©tr(T> 2 ) ©tr(_D 3 ) = 
0 = 1, and hence conclude that the conditions of Theorem |7| are fulfilled. 
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Furthermore, we successively obtain 


D* 


I ®D®D 2 = 

(dT AD*)~ = 


0 -2 1 \ 

2 0 3, AD* 

-1 - 30 / 

°\ 

2 1 , D*(dAD*Y 


AD*(d~AD*)- 



4 2 5 \ 

5 3 6, 
1 - 12 / 



The minimum value of the objective function, which shows the minimal 
violation of the due dates, is given by A = ((AD*(d~AD*)~)~ d) 1 / 2 = 2. 

The optimal schedule has the latest start and finish times defined by the 
vectors 


a; = A D*(d-AD*)~ 



y = A AD*(dAD*)~ 



5.2 Minimization of maximum deviation of finish times 

After rewriting problem Q in terms of the semifield lR m ax,+, the problem 
becomes 

n n 

minimize 

i =1 j =1 

n n 

subject to (J) a ij x 3 = 2/i> © bijXj < x t , 

3 = 1 i =1 

Vi<fi , i = l,...,n. 


© c l3 y 3 < a*, 
i=i 


In addition to the previously introduced matrix-vector notation, we de¬ 
fine the vector f = (fi) to write the problem 


rp _ 

minimize 1 yy 1, 

subject to Ax = y, Bx < x, Cy < x, y < /. 


( 12 ) 


The following result provides a solution to the problem. 

Theorem 8. Let A be a regular matrix, the matrix D = B © CA satisfy 
the condition Tr (D) < 1, and f be a regular vector. Then, the minimum 
value in problem m is equal to 

A = (AD*(1 t AD*)~)-1, 
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and attained if 

x = uD*(l T AD*)~, y = aAD*(l T AD*)~, a < (f~ AD* (\ T AD*)~)~ l . 


Proof. As before, we substitute y = Ax and combine the first two inequal¬ 
ities into one inequality Dx < x, where D = B ® CA. This inequality is 
then solved by using Theorem [2] to give the result x = D*u, where u is a 
regular vector. 

Furthermore, we write the last inequality constraint as AD*u < /, and 
apply Lemma |T] to find the solution u < ( f~AD *)~The problem takes 
the form 

minimize 1 T AD*u(AD*u)~l, 
subject to u < (f~AD*)~. 

First, we remove the constraints and solve the obtained unconstrained 
problem. By applying Theorem @1 where both matrices A and B are re¬ 
placed by AD* , and both vectors p and q by 1, we find the minimum A = 
(AD*(1 T AD*)-)~1, which is attained at the vector u = a(l 1 AD*)~ , 
where a > 0. 

To find the values of the parameter a, which meet the condition u < 
(f~AD*)~ , we solve the inequality a(l 7 AD*)~ < (f~AD*)~ . By apply¬ 
ing Lemma[T]with ot as the unknown, we have a < (f~AD*(1 J AD*)~)~ l . 

It remains to turn back to the vectors x and y , and thus complete the 
proof. □ 

Example 2. Suppose that we need to minimize the maximum deviation of 
finish times in the project from Example [TJ, where, instead of the due dates, 
deadline constraints apply, given by the vector 



We take advantage of intermediate results of the previous example to 
obtain 


-5 

(l r AD*)~ = I -3 


-5 


D*(1 t AD*)~ =| -3 I , 


-6 

AD*(1 t AD*)~ = 


-6 


-1 

0 

-4 


According to Theorem [81 the minimum deviation of finish times, which 
can be achieved in the project, is equal to A = (AZD*(1 7 AD*)~)~1 = 4. 
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The optimal schedule is provided by the vectors 


x = aD*(l T AD*)~ = a 



y = aAD*(l T AD*)~ = a 


-1 

0 

-4 


where the condition a < (f-AD*{l T AD*)-)~ 1 = 6 must be satisfied. 

In terms of standard operations, the elements of the vectors x = (x\, X 2 , x%) T 
and y = (yi, 2 / 2 , 2 / 3 ) T are written as 


x\ = a — 5, X 2 = ol — 3, £3 = a — 6 , 

yi = a- 1, 2/2 = «, y 3 = a-4, a < 6 . 


5.3 Minimization of maximum flow-time 

Consider problem (]3|), and rewrite it in terms of the semifield lR m ax,+ - As a 
result, we obtain the problem 

n 

minimize e Xi 1 Vu 
2=1 

n n n 

subject to (J) dijXj = y i: ^ bijXj < x i} ^ c^y^ < x t , 
j =1 3 =1 i=i 

9 i<Xi, i = 1,... ,n. 

Furthermore, we add the vector g = {gt) ■ Switching to matrix-vector 
notation puts the problem in the form 

minimize x~y, 

y (13) 

subject to Ax = y , Bx < x, Cy < x, g < x. 

A complete solution of the problem is given as follows. 

Theorem 9. Let A be a matrix with spectral radius A > 0, and the matrix 
D = B © CA satisfy the condition Tr(D) < 1. Then, the minimum value 
in problem m is equal to 

n —1 

e = A © 0 0 ti J / fc (AD h ■ ■ ■ AD ik ), 

k =1 1<hH- \-ik<n—k 

and all regular solutions are given by 

x = {0~ l A 0 D)*u, y = A{0~ 1 A© D)*u, u>g. 
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Proof. By substitution of the equality constraint y = Ax , we eliminate the 
vector y. Then, we combine all inequality constraints into one to write the 
problem 

minimize x~Ax, 

subject to (B ® CA)x © g < x. 

This problem has the form of that at Q, where B is replaced by D = 
B © CA. Thus, a direct application of Theorem [5] yields the required 
solution. □ 


As an illustration of the solution obtained, we present the next example. 

Example 3. Consider the problem of minimizing the maximum flow-time in 
the project, which has the start-finish, start-start and finish-start constraints 
defined as in Example [H and release time constraints given by 



First, we verify the conditions of Theorem [9j We calculate the matrices 


A 2 = 




A 3 


12 8 3 \ 

9 9 5, 

8 4 6 / 


and then evaluate the traces to find A = tr A © tr 1 / 2 (A 2 ) © tr 1 / 3 (A 3 ) = 4. 
Since A > 0 = 1 and Tr(Z?) = 0 = 1, the conditions of the theorem are 
satisfied. 

To evaluate the minimum 0 of the objective function, we calculate the 
matrices 


AD 


/ 3 

2 

5 \ 

/ 4 

1 

5 

4 

2 

6 ’ 

AD 2 = 5 

2 

5 

V 1 

-2 

1 / 

\ o 

-1 

2 


ADA 


7 5 7 \ 

8 5 8, 
5 13/ 


A 2 D = 


7 6 9 \ 
7 5 9 
3 2 5 / 


After taking the traces and considering A, we have the result of minimiz¬ 
ing the maximum flow-time, given by 9 = Affitr(AZAffi AZA 2 )©tr 1 / 2 (AZAA© 
A 2 D) = 4. 

To describe the solution offered by Theorem [9j we need the matrices 

/ 0 -2 1 \ / 0 -2 1 \ 

9- 1 A®D=l 1-1 3 , (0~ 1 A © D) 2 = I 2 -1 2 . 

V —1 —6 —2 / \ -1 -3 0 / 
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By combining these matrices with the identity matrix, we obtain 


(0~ A ® D)* = 


0 -2 1 
2 0 3 

-1 -3 0 


Note that the last matrix can be represented in the form 
(O^A^D)* = [ 3 ] ( -1 -3 0 ) 

The vector of optimal start times provided by Theorem [U] is given by 
x = (9~ V A © D)*u = | 3 J ( — 1 — 3 0 ) u, u> g. 


To simplify the solution, we introduce the new variable v = ( —1 —3 0 ) u, 

which has to satisfy the condition v = ( — 1 —3 0 ) u > ( — 1 —3 0 ) g = 

1. 

The vectors of optimal start and finish times now become 



Using standard operations yields the elements of the vectors given by 

X\ = v + 1 , x 2 = v + 3, x 3 = v, 

y\ = V + 5, y 2 =v + 6, y 3 = v + 2, V > 1. 


5.4 Minimization of makespan 

In the framework of the semifield !R m ax,+ > problem (J4]) of minimizing the 
makespan is rewritten as 

n n 

minimize ©2/i©^- 1 , 
i= 1 3 =1 

n 

subject to (J) aijXj = y h gi <x t < hi, 
j =i 

Vi<fi, i = l,...,n. 

By adding the vector h = (hi) and using 1 to indicate the vector of 
ones, we represent the problem as follows: 


rp _ 

minimize 1 yx 1, 

subject to Ax = y, g < x < h, y < f. 


(14) 
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Theorem 10. Let A be a column-regular matrix, h and f be regular vec¬ 
tors such that {h~ © f A)g < 1. Then, the minimum value in problem 
dD is equal to 

6 = 1 T A(I © gh~)l, 
and all regular solutions are given by 

x = (J © A)u , y = A(I © 6~ 1 11 T A)u, 

where 

g < u < (( h~ © f~A)(I © Q^ll T A))~. 

Proof. As before, we first substitute y = Ax. Solving the inequality Ax < 
f by Lemma [T] yields x < ( f~A)~. Then, we take the upper boundaries 
x < h and x < ( f~A )~, and apply conjugate transposition to rewrite the 
inequalities as x~ > h~ and x~ > f A. By coupling both inequalities 
into one, and again taking the conjugate transposition, we obtain one upper 
bound x < (h~ © f~A)~ . 

We write the objective function as l 1 Axx l = x~ll T Ax to obtain 
the problem 

minimize x 

(15) 

subject to g < x < (h © / A) . 

The problem obtained is of the form of (1101) . where A is replaced by 
11 t A and h by (h~ © f~A)~. To apply Theorem [HI we need to find 
the spectral radius of the matrix 11^ A. We first calculate, for each k = 

(ll T A) fc = (l T Al) fe_1 ll T A, tr(ll T A) fc = (l T Al) fc , 

from which it follows that the spectral radius is equal to A = 1 T A1 > 0. 
Furthermore, we consider the minimum given by 

n— 1 

6 = 1 t A1 © (1 T A1) 0((l T Al) _1 h _ ll r Ag) 1/fc . 

fc=l 

Suppose that l 1 Al < h~ll 7 Ag. Since the inequality (l 1 Al) -1 h~ll 7 Ag > 
1 holds, we have ((1 T Al) _1 h _ ll T Agi) 1 / fc < (1 T Al)- 1 h~ll 1 Ag, and, 
therefore, conclude that 0 = h~ll 7 Ag. 

On the other hand, if l 7 Al > h ll 1 Ag , then we have 6 = l 7 Al. 

By combining both results, we finally obtain 9 = l 7 Al © h~ll T Ag = 
l T A{I®gh~)l. 

To describe the solution set according to Theorem 0 we examine the 
matrix 

n—l n— 1 

(0 -1 11 t A)* = 0(0" 1 ll T A) fc = I © 0- 1 0(0- 1 l T Al) fc - 1 ll r A. 

k =0 k =1 
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Considering the inequality 9 > l 7 Al, we have (0^ 1 ll r A)* = 7® 
9~ 1 11 I A. Substitution into the solution provided by Theorem [6] yields 
x = (7 © 0 -1 11 t A) ii, where the vector u satisfies the condition g < u < 
((h~ © f~A)(I © 0~ l ll T A))~ . 

Finally, we represent the vector y = Ax, which completes the proof. □ 

Example 4. Assume that we need to find a schedule with the minimum 
makespan under the start-finish, release time and deadline constraints, which 
are defined by the matrix A and the vectors g and f in the previous 
examples. Suppose that, in addition, release deadlines have to be taken into 
account, given by the vector 



To verify that the condition of Theorem 1101 is fulfilled, we calculate the 
vectors 


ra=(- 2 -3 -4), h-®f-A=(- 2 -3 -2). 

and then obtain the required condition in the form ( h~ © f A)g = 0 = 1. 

We now apply Theorem [10] to find the optimal schedule. The calculation 
of the minimum makespan 9 involves 

( 0 -1 0 \ ( 0 \ 

I@gh~ = -1 0 0, (I(Bgh~)l = I 0 I , 1 T A = ( 4 3 2 ), 

1-2-20/ W 

which results in 9 = l 7 A(7 © gh~) 1 = 4. 

To represent the solution, we need to find the matrices 

/ 4 3 2 \ / 0 —1 — 2 \ 

11 t A= 4 3 2, 7 © 9~ 1 11 T A = I 0 0 -2 , 

\ 4 3 2 / \ 0 -1 0 / 

and to form the vectors 

«i = 9 = ^ 2 j , u 2 = ((h~ © f~A)(I © 9- 1 11 t A))~ = ^ 3 j . 


The vectors of the optimal start and finish times in the project take the 
form 


-in T 


/ 0 

-1 

- 2 \ 

( 4 

3 

2 \ 

0 

0 

-2 u , 

y = Ax = 3 

3 

1 

\ 0 

-1 

0 ) 

V 2 

1 

2 / 
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where the vector u satisfies the condition u\ < u < m 2 ■ 

Note that the solution can be simplified as follows. The lower and upper 
bounds Mi and M 2 for the vector u yield the corresponding bounds on x 
in the form 


®i = (I © 6 l ll 7 A)u\ = ^ 2 ^ , x 2 = (I @0 1 11 T A)m 2 = ^ 3 ^ . 

Since the first and third elements of X\ and X 2 coincide, the elements 
of the vector x can be directly defined by aq = 2, 2 < :r 2 < 3 and X 3 = 2. 

After calculating the vector y = Ax , we write the results using a pa¬ 
rameter v and standard operations as follows: 

x\ = 2, X 2 = v, x 3 = 2 , 

2/1 = 6, 2/2 = v + 3, 2/3 = 4, 2 < v < 3. 

6 Conclusions 

In this paper, we demonstrated new solution techniques, based on the mod¬ 
els and methods of tropical mathematics, for a class of project scheduling 
problems with minimax objectives. The paper extends and generalizes re¬ 
sults in [21, 2U [271 26, [29] by solving problems not previously considered, 
which involve different objective functions and/or more complicated systems 
of constraints. 

It was shown that many typical constraints and objectives that occur 
in time-constrained project scheduling are naturally represented in terms 
of tropical algebra in a compact vector form. Specifically, the start-finish, 
start-start and finish-start precedence relationships are readily given by lin¬ 
ear vector equations and inequalities, whereas the project makespan, the 
maximum deviation from due dates and the maximum flow-time optimiza¬ 
tion criteria can be well written as nonlinear functions defined on vectors 
through a multiplicative conjugate transposition operator. 

The scheduling problems of interest were reduced to tropical optimiza¬ 
tion problems, which are solved through the application and further devel¬ 
opment of recent results in tropical optimization. The solutions are given in 
direct, explicit vector forms ready for formal analysis and straightforward 
computation with low polynomial complexity. The explicit form of the so¬ 
lutions can be considered as a definite advantage over the algorithmic tech¬ 
niques normally used in scheduling. Specifically, the problems under study 
can be formulated as linear programs, which offers algorithmic solutions us¬ 
ing one of the iterative computational schemes of linear programming, but 
does not guarantee a direct closed-form solution. 
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The contribution of the paper is twofold. First, we have developed new 
applications of tropical optimization, formulated new optimization prob¬ 
lems, extended existing solution methods to these problems, and derived 
direct solutions. Second, we have devised a new solution approach for time- 
constrained scheduling problems, which involves the representation of the 
problems in the tropical mathematics setting, and the application of meth¬ 
ods from tropical optimization to obtain direct solutions in a closed form. 
Analytical techniques proposed and described in the paper can serve as a 
template for the solution of other optimization problems, and for the appli¬ 
cation of the solutions to real-world problems in various fields. 

Possible lines of further investigation include the extension of the ap¬ 
proach to account for new types of constraints and objectives in the schedul¬ 
ing problems under study, and to solve new classes of problems, including 
scheduling problems with renewable resources. Computational experiments 
with real-world data present another research topic of interests. 


References 

[1] M. Akian, R. Bapat, and S. Gaubert. Max-plus algebra. In L. Hog- 
ben, editor, Handbook of Linear Algebra , Discrete Mathematics and Its 
Applications, pages 25-1-25-17. Taylor and Francis, Boca Raton, FL, 
2007. doi:10.1201/9781420010572.ch25, 

[2] A. Aminu and P. Butkovic. Non-linear programs with max-linear con¬ 
straints: A heuristic approach. IMA J. Manag. Math., 23(l):41-66, 
2012. doi:10.1093/imaman/dpq020. 

[3] C. Artigues and F. Roubellat. A polynomial activity insertion algo¬ 
rithm in a multi-resource schedule with cumulative constraints and 
multiple modes. European J. Oper. Res., 127(2):297—316, 2000. 
doi:10.1016/S0377-2217(99)00496-8, 

[4] F. L. Baccelli, G. Cohen, G. J. Olsder, and J.-P. Quadrat. Synchroniza¬ 
tion and Linearity. Wiley Series in Probability and Statistics. Wiley, 
Chichester, 1993. 

[5] J.-L. Bouquard, C. Lente, and J.-C. Billaut. Application 
of an optimization problem in max-plus algebra to schedul¬ 
ing problems. Discrete Appl. Math., 154(15):2064-2079, 2006. 
doi:10.1016/j.dam.2005.04.Oil. 

[6] P. Butkovic. Max-linear Systems. Springer Monographs in Mathemat¬ 
ics. Springer, London, 2010. doi : 10.1007/978-1-84996-299-5. 

[7] P. Butkovic and K. P. Tam. On some properties of the image set 
of a max-linear mapping. In G. L. Litvinov and S. N. Sergeev, 


24 



editors, Tropical and Idempotent Mathematics, volume 495 of Con- 
temp. Math., pages 115-126. American Mathematical Society, 2009. 
doi:10.1090/conm/495/09694. 

[8] R. Cuninghame-Green. Minimax Algebra, volume 166 of Lecture Notes 
in Economics and Mathematical Systems. Springer, Berlin, 1979. 
doi:10.1007/978-3-642-48708-8. 

[9] R. A. Cuninghame-Green. Describing industrial processes with inter¬ 
ference and approximating their steady-state behaviour. Oper. Res. 
Quart., 13(1):95—100, 1962. doi : 10.2307/3007584, 

[10] B. De Reyck and W. Herroelen. A branch-and-bound procedure for 
the resource-constrained project scheduling problem with generalized 
precedence relations. European J. Oper. Res., 111(1): 152—174, 1998. 
doi:10.1016/S0377-2217(97)00305-6, 

[11] E. L. Demeulemeester and W. S. Herroelen. Project Scheduling, vol¬ 
ume 49 of International Series in Operations Research and Management 
Science. Springer, New York, 2002. doi : 10.1007/bl01924. 

[12] M. Fiedler, J. Nedoma, J. Rann'k, J. Rohn, and K. Zimmermann. Lin¬ 
ear Optimization Problems with Inexact Data. Springer, Berlin, 2006. 
doi:10.1007/0-387-32698-7. 

[13] B. Giffler. Scheduling general production systems using sched¬ 

ule algebra. Naval Res. Logist. Quart., 10(1):237—255, 1963. 
doi:10.1002/nav.3800100119. 

[14] J. S. Golan. Semirings and Affine Equations Over Them , volume 556 
of Mathematics and Its Applications. Kluwer Acad. Publ., Dordrecht, 
2003. doi:10.1007/978-94-017-0383-3. 

[15] M. Gondran and M. Minoux. Graphs, Dioids and Semirings, volume 41 
of Operations Research/ Computer Science Interfaces. Springer, New 
York, 2008. doi:10.1007/978-0-387-75450-5. 

[16] B. Heidergott, G. J. Olsder, and J. van der Woude. Max Plus at 
Work. Princeton Series in Applied Mathematics. Princeton Univ. Press, 
Princeton, NJ, 2006. 

[17] A. J. Hoffman. On abstract dual linear programs. Naval Res. Logist. 
Quart., 10(l):369-373, 1963. doi : 10.1002/nav. 3800100131. 

[18] L. Houssin. Cyclic jobshop problem and (max,plus) algebra. In S. Bit- 
tanti, A. Cenedese, and S. Zampieri, editors, Proceedings of the 18th 
IFAC World Congress, 2011, volume 18 of World Congress, pages 2717- 
2721. IFAC, 2011. doi:10.3182/20110828-6-IT-1002.03095. 


25 



[19] V. N. Kolokoltsov and V. P. Maslov. Idempotent Analysis and Its Ap¬ 
plications, volume 401 of Mathematics and Its Applications. Kluwer 
Acad. Publ., Dordrecht, 1997. doi : 10.1007/978-94-015-8901-7. 

[20] A. A. Korbut. Extremal spaces. Soviet Math. Dokl., 6(5):1358-1361, 
1965. 

[21] N. Krivulin. Algebraic solutions to scheduling problems in project man¬ 
agement. In I. J. Rudas, A. Zaharim, K. Sopian, and J. Strouhal, ed¬ 
itors, Recent Researches in Communications, Electronics, Signal Pro¬ 
cessing and Automatic Control, pages 161-166. WSEAS Press, Athens, 
2012. 

[22] N. Krivulin. A solution of a tropical linear vector equation. In S. Yenuri, 
editor, Advances in Computer Science, volume 5 of Recent Advances in 
Computer Engineering Series, pages 244-249. WSEAS Press, Athens, 

2012. arXiv:1212.6107, 

[23] N. Krivulin. Explicit solution of a tropical optimization problem with 
application to project scheduling. In D. Biolek, H. Walter, I. Utu, 
and C. von Lucken, editors, Mathematical Methods and Optimization 
Techniques in Engineering, pages 39-45. WSEAS Press, Athens, 2013. 
arXiv:1303.5457, 

[24] N. Krivulin. A constrained tropical optimization problem: Complete 
solution and application example. In G. L. Litvinov and S. N. Sergeev, 
editors, Tropical and Idempotent Mathematics and Applications, volume 
616 of Contemp. Math., pages 163-177. AMS, Providence, RI, 2014. 
arXiv:1305.1454, doi:10.1090/conm/616/12308. 

[25] N. Krivulin. Tropical optimization problems. In L. A. Petrosyan, J. V. 
Romanovsky, and D. W. K. Yeung, editors, Advances in Economics 
and Optimization, Economic Issues, Problems and Perspectives, pages 
195-214. Nova Sci. Publ., New York, 2014. arXiv: 1408.0313. 

[26] N. Krivulin. Extremal properties of tropical eigenvalues and solutions 
to tropical optimization problems. Linear Algebra Appl, 468:211-232, 
2015. arXiv:1311.0442, doi:10.1016/j.laa.2014.06.044. 

[27] N. Krivulin. A multidimensional tropical optimization prob¬ 

lem with nonlinear objective function and linear constraints. 
Optimization, 64(5): 1107—1129, 2015. arXiv: 1303.0542, 

doi:10.1080/02331934.2013.840624. 

[28] N. Krivulin. Tropical optimization problems in project scheduling. In 
Z. Hanzalek, G. Kendall, B. McCollum, and P. Sucha, editors, MISTA 
2015 Proceedings, pages 492-506. MISTA, 2015. arXiv: 1502.06222. 


26 



[29] N. Krivulin. Tropical optimization problems with application to project 
scheduling with minimum makespan. Ann. Oper. Res., pages 1-18, 
2015. doi:10.1007/s 10479-015-1939-9. 

[30] N. K. Krivulin. Solution of generalized linear vector equations in idem- 
potent algebra. Vestnik St. Petersburg Univ. Math., 39(l):16-26, 2006. 

[31] G. Litvinov. Maslov dequantization, idempotent and tropical mathe¬ 
matics: A brief introduction. J. Math. Sci. (NY), 140(3) :426-444, 2007. 
arXiv:math/0507014, doi:10.1007/sl0958-007-0450-5. 

[32] K. Neumann, C. Schwindt, and J. Zimmermann. Project Scheduling 
with Time Windows and Scarce Resources. Springer, Berlin, 2 edition, 
2003. doi:10.1007/978-3-540-24800-2. 

[33] S. N. N. Pandit. A new matrix calculus. J. SIAM, 9(4):632—639, 1961. 
doi:10.1137/0109052, 

[34] A. A. B. Pritsker, L. J. Watters, and P. M. Wolfe. Multiproject schedul¬ 
ing with limited resources: a zero-one programming approach. Man¬ 
agement Sci., 16(1):93-108, 1969. doi : 10.1287/mnsc . 16.1.93, 

[35] I. V. Romanovskii. Asymptotic behavior of dynamic programming pro¬ 
cesses with a continuous set of states. Soviet Math. Dokl., 5(6): 1684 
1687, 1964. 

[36] D. Speyer and B. Sturmfels. Tropical mathematics. Math. Mag., 
82(3):163-173, 2009. 

[37] V. T’kindt and J.-C. Billaut. Multicriteria Scheduling. Springer, Berlin, 
2 edition, 2006. doi : 10.1007/bl06275, 

[38] M. Vanhoucke. Project Management with Dynamic Scheduling. 
Springer, Berlin, 2012. doi : 10.1007/978-3-642-40438-2, 

[39] N. N. Vorob’ev. The extremal matrix algebra. Soviet Math. Dokl., 
4(5):1220-1223, 1963. 

[40] K. Zimmermann. Some optimization problems with extremal opera¬ 
tions. In B. Korte and K. Ritter, editors, Mathematical Programming 
at Oberwolfach II, volume 22 of Mathematical Programming Studies, 
pages 237-251. Springer, Berlin, 1984. doi : 10.1007/BFb0121020. 

[41] K. Zimmermann. Disjunctive optimization, max-separable problems 
and extremal algebras. Theoret. Comput. Sci., 293(l):45-54, 2003. 
doi:10.1016/S0304-3975(02)00231-1, 


27 


[42] U. Zimmermann. Linear and Combinatorial Optimization in Ordered 
Algebraic Structures, volume 10 of Annals of Discrete Mathematics. 
Elsevier, Amsterdam, 1981. 


28 



